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0\ \ Abstract: We study the compatibility of a light bottom squark M^ < 0{lQGeV) 

(^ ! in the unconstrained MSSM. We consider the one-loop radiative corrections which 

O 



are large for a heavy gluino {>0{lhQGeV)). We then consider the renormalization 

group flow up to the Grand Unified scale. For most regions of the parameter space 

, with a light sbottom we find colour and charge breaking minima. Only a small 

D , window in gluino mass and tan (5 is consistent with all bounds. This is alleviated by 



a light gluino, which is however only marginally experimentally allowed. 
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1. Introduction 

The experimental bound from LEP on the lightest supersymmetric particle (LSP), 
assuming it is a neutralino, is given by 0] 

M^o>40.9GeV; (OPAL), (LI) 

and similar numbers from DELPHI {36.7 GeV) [|, L3 (38.2 GeX^) || and ALEPH 
[37.5 GeV) [Q. This bound assumes the supersymmetric grand unified relation be- 
tween the gaugino masses at the weak scale: Mi = (5/3) tan^ 9w M2 (where 9w is 
the electroweak mixing angle) and employs the chargino search. In a recent paper 
[^] it was shown that if you drop this theoretical assumption a LSP neutralino even 
as light as 34MeV is consistent with all experiments.^ 

It is the purpose of this letter to do a similar study for light bottom squarks. 
Light top squarks have been extensively studied elsewhere 0]. For large values of 
tan P (the ratio of the vacuum expectation values of the two neutral CP-even Higgs 
bosons in the MSSM) it can be natural to have light bottom squarks as well, as we 
discuss in more detail below. Both the DO and CDF experiments have performed 
direct searches for the lightest bottom squark [|, |] obtaining 

m^>115GeV, (DO) (L2) 

m^>146GeV, (CDF). (L3) 



^For a discussion of astrophysical bounds see n S. 



We have given the maximum bound which is obtained for a vanishing neutrahno LSP 
mass. In general the bound depends on the LSP mass and becomes less sensitive 
as the mass difference between the LSP and the squark is decreased. For smaller 



mass differences and also for smaller bottom squark masses the LEP searches [|I0, 



Tl| , [T^, |T^, ^ are more sensitive. However, even in this case there remains a gap 
at very small mass differences to the LSP which becomes more pronounced for very 
small squark masses m^ < 0(10 GeV) [jl4|. At such low masses the decay b -^ bxi is 
kinematically suppressed by the final state quark even for vanishing neutralino mass. 
A dedicated search for the top squark with a small mass difference (AM) to the LSP 
has been performed [^ reaching as low as AM = l.QGeV; the threshold for the 



decay t -^ x°c. We are not aware of such a search for light bottom squarks. 

Light squarks can directly contribute to the hadronic cross section at the Z° peak. 
As we will discuss below, the doublet and singlet squarks mix and for specific mixing 
parameters the coupling to the Z^ can even vanish. Thus this constraint restricts 
the range of sbottom mixing but can not exclude a light sbottom. This constraint 
turns out to be very mild since the light sbottom is dominantly an SU{2)l singlet. 

Very light bottom squarks have recently been investigated in Refs. [^, jl^. In 
r^ a possible influence on the parameter R{s) = a{e^e~ -^ hadrons)/(j(e"'"e~ — > 



/i"'"/i~) was studied. For a 6-squark the asymptotic contribution is only 1/12. This is 
1/4 that of a b-quark due to the missing spin degeneracy and is below the experimen- 
tal sensitivity [|1^. In [^ the effect of a light b-squark on the electroweak precision 



data and on the MSSM Higgs sector was investigated. It was found to be consistent 
with the precision data provided the scalar top quark is not too heavy. The upper 
bound on the lightest CP-even Higgs boson in the MSSM is slightly lowered. 

In the following we discuss the theoretical implications of a very light bottom 
squark. We focus on the embedding into the MSSM. We first study the pole mass of 
the bottom squark at one-loop, including in particular radiative corrections from the 
gluino which are large, and also corrections from top and bottom Yukawa couplings 
to the Higgs boson masses. We then study the renormalization group flow of the 
right-handed bottom squark mass squared for both universal and non-universal scalar 
fields at the GUT scale. We consider the constraints from colour and charge breaking 
minima (CCB). The constraints are relaxed by a light gluino. We finish with a brief 
discussion of bounds on a light gluino, before we conclude. 

2. Parameters and Constraints 

In the MSSM there are two bottom squarks. The SU{2)l current eigenstates are 
denoted 6^, Br, where bi is a doublet squark and 6r is a singlet squark. The corre- 
sponding states for the top squark are t^, t^. The mass matrix of these squarks in 
the current eigenstate basis is given for example in [1191 ^^'^ ^^^ one-loop radiative 



corrections are given in [20 . 



The mass eigenstates depend on the following parameters in the standard MSSM 



notation [^]: M^^ j^, the doublet and singlet soft-breaking squark masses, respec- 



tively; Mw,z, the gauge boson masses; tan/3, the ratio of the vacuum expectation 
values of the two neutral CP-even Higgs fields; A^, At, the tri-linear soft breaking 
terms; /i, the Higgs mixing parameter; and Mg, the gluino mass. In the following rrib 
and nit denote the bottom and top quark mass, respectively. We shall denote the 



lighter bottom squark 62 and the heavier one 61 in accordance with pO[. The scalar 
bottom mixing angle we denote O^. All the above parameters are considered to be 
DR-Tunning parameters. 

Besides the direct searches we have discussed in the introduction a light b-squark 
would also contribute to the hadronic cross section at the Z° peak, 0"°^^. The exper- 



imental bound for any contribution beyond the SM is [18, p2| 



2 



Aal,{Z')<0.U2nb, (2a). (2.1) 

At tree-level this requires the sbottom mixing angle to lie in the range (for sin 6'^ = 
0.2315 and N^ = 3) 

I sin ^^1 < 0.535. (2.2) 

At 2a, zero mixing is consistent with the data. In our analysis we have included the 
one-loop contribution to 0"°^^^ from the scalar bottom. We have only plotted points 
which are consistent with the bound ( |2.2| ). It turns out that this constraint has no 
effect on Figs. |I|-0. For 62 satisfying ( |2.2| ) the heavier bottom squark, 61, couples 



unsuppressed to both the photon and the Z^. In order to avoid experimental bounds 
from LEPl and LEP2 we must therefore require 



mi^^200GeV. (2.3) 

In our scans below, we shall employ this bound, as well. 

A light sbottom contributes to the running of the strong coupling a^ between rrir 
and Mzo. In order to see whether this is consistent with the data one must include 
the light sbottom both in the determination of a^ in a given experiment and also 
in the beta function. This is beyond the scope of this letter. However, it has been 
performed for a light gluino [^, ^. The most recent study [0] with the smallest 
experimental error in a^ was able to exclude a light gluino at the 70% C.L.. The 
contribution to the beta function at one-loop of a singlet sbottom is 1/12 that of 
a gluino. We thus expect the effect to be significantly smaller and beyond present 
experimental sensitivity. 



^We note that the Standard Model prediction of cr'^^^ is currently 1.7 a below the measured 
value. In principle this could be exactly compensated by the light bottom squark p^ . In the 
following we choose to focus only on the experimental upper bound on a new contribution. 
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Figure 1: a) Contour plot of the tree level light sbottom mass mj^ as a function of the 
right handed singlet soft SUSY breaking mass M^ and tan/3. The tree level mass-squared 
sign^irij^ )|r?T,? \^'^ is getting negative for tan/3 > 10 and for various values of M^ indicated 
in the figure. The small dashed band indicates values of the light sbottom in the region 
1 — 10 GeV. The soft trilinear coupling Af, has been set to zero in this plot, b) As in Fig. la 
but for the physical 1-loop mass. Sbottom masses of 50, 75, 100, 125 GeV are indicated 
for comparison. In the large shaded region, the running DR mass squared (M^ ) is getting 
negative at a scale which lies below the GUT scale. 



3. Sbottom Pole Mass 



We now investigate the effect of radiative corrections on the sbottom pole mass. 
mi^ depends at tree-level on the parameters Mq, Mfj, Af,, fi and tan/3. At one- 
loop ||20|, there is a further dependence on the stop mass parameters, Mg, and Ma^ 
The dependence on the stop sector and the Higgs sector parameters is weak and 
we fix them to At = 300 GeV, M^^ = 300 GeV and Ma = 400 GeV. We also fix 
the following SM parameters at the Z^ scale to: mt{pole) = 175 GeV, sin^ 6^ = 
0.2315, nib = 2.9 GeV, and a^ = 0.12. As we discuss now, the dependence on M^ , 
Ab and on Mg is strong. 

The first case we examine is that of a heavy gluino of 200 Gel^, just above the 
current experimental bound of 180 GeV [^. For now, we fix the remaining input 
parameters: M^^ = 250 GeV, A^ = GeV, (x = 250 GeV, at the Z° scale. In Fig. |I|a,b 
we present contour plots of the lightest bottom squark mass, m^^ in the (M^ -tan/3) 
plane. We display both tree level and physical 1-loop pole masses. The narrow 



^We do not include chargino-neutralino corrections since they are small [g^. Also, the light 
Higgs mass has been set to Mz- Variation of the Higgs mass to its upper limit affects very weakly 
the light sbottom mass. 
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Figure 2: a) Contour plot of the tree level light sbottom mass as a function of the right 
handed singlet soft SUSY breaking mass Mj^ and the trilinear coupling Ai,{Mz)- The 
value of tan/3 is fixed to 15. 2b) The same for the physical 1-loop light sbottom pole mass. 
The light sbottom mass contours, 1 — 10 GeV of this figure is completely within the shaded 
region of Fig. lb. 



shaded strip corresponds to masses in the range 1 — 10 GeV . The area to the left of 
this narrow strip in Fig. lb is excluded since the scalar bottom pole mass squared 
turns out to be negative. The effect of the radiative corrections is significant for 
tan/3 ~ 15. One can see up to a 40 GeV difference between the tree level and the 
1-loop physical mass. They tend to push a fixed sbottom mass to larger values of 
Ml . The solution region for a light sbottom is quite narrow and somewhat fine- 
tuned. A variation of 1 GeV of M^ results in a variation of more than 5 GeV in 
the light sbottom mass in the 0{< lOGeV^) region. Thus when determining the 
supersymmetric parameters for a light bottom squark the radiative corrections need 
to be taken into account. 

With the above input values we obtain for the other (physical 1-loop) masses : 



m^^ = 255 - 300 GeV, m^^ 



191 - 228 GeV, mi^ = 394 - 416 GeV. All these masses 



satisfy the current experimental bounds. 

As we have already mentioned above, another parameter which plays a crucial 
role in determining the mass of the bottom squark is the trilinear coupling Ah{Mz)- 
This parameter enters in both tree level and 1-loop sbottom mass corrections and 
the effect of its variation is presented in Fig.^. The input parameter tan/? is fixed 
to tan/5 = 15. As we can see, the one loop radiative corrections shift the mass 
contours by at most tens of GeV. Large values of Aiy{Mz) typically give smaller 
sbottom masses. The other physical masses (m^ , m^^, rrif:^) vary inside the region 
we mentioned in the previous paragraph. 
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Figure 3: a) The same as in Fig. lb 3b) and in Fig. 2b but for a light gluino of mass, 
Mg = 3 GeV. The tree level results are those given in Fig. la and Fig. 2a . 



Further important radiative correction to the bottom squark masses are those 



arising from the loops involving gluinos |]20|. This is obvious when one compares the 
Figs. 1^, where the gluino mass is taken to be 200 GeV with Fig. ^ where its mass is 
set to be 3 GeV . We see that in the (M^ -tan /?) plane there is almost no effect from 
radiative corrections for a light gluino compared to the tree-level result presented in 
Fig.|l|a. In the (M^ -A\^ plane the effect is much less dramatic for a light gluino, 
although there is still a qualitative difference in the A^, dependence compared to the 
tree-level result. 

In summary, a light bottom squark can be consistently implemented at one- 
loop in the MSSM. The effect of radiative corrections as a function of Mg and A;, 
is substantial, up to several tens of GeV and must be considered when determining 
the supersymmetric parameters. 



4. Renormalization Group 



4.1 Universal Scalar Masses at Mx 

We next consider the embedding of the MSSM in a more unified theory at a high 
scale, Mx = 0{1Q^^ GeV). Having extracted the DR quantity M^ from the Figs |I], 
^,|^ we would thus like to see if these values are compatible with the renormalization 
group running of the MSSM up to Mx- In order to qualitatively understand the 
evolution we first discuss an approximate analytic solution for the M^ running mass. 
We present the full numerical analysis below. The renormalization group evolution 



of the M^ mass is given |2^ by, 



where 



IQtt 



2 feff 

dt 



^y,'n - f ^|M| - ^^glMl + IgfTriYm^) 



(4.1) 



nf 



Tr(Fm2) ^ M^^ - M^^ + 5:(M?^ - 2M,^ + M. - M?^ + M,|^ 



i=l 



(4.2) 
(4.3) 



To start, we assume that the contribution from the bottom Yukawa couphng is small, 
i.e. we restrict ourselves to the region tan/3 ~ 10. We also assume all of the squark, 
slepton and Higgs-boson masses are the same at the GUT scale, i.e. TiiYim?) remains 
zero at all scales. In the case where the gluino is heavy, Mg = 200 GeV, the running 



of M^^ "freezes" below Mg [g|] to 



M 



bR 



ml + C3 + -Ci - - sin^ e^Ml cos(2/5) 



(4.4) 



where mo is the common squark, slepton and Higgs-boson mass at the GUT scale 
and 



Cif,) = --Mi 



1- 



a?(M 



X 



Csif^) 



9" 



-M' 



af(/i) 

«i(/i) 



(4.5) 



The D-term contribution, — ^ sin^ O^M"^ cos(2/?) in Eq. (^^) , is positive since cos(2/5) 
is negative, and also the bino contribution Ci is positive. The dominant term in 
the evolution of M? is the gluino contribution, C3. Now, suppose that all the scalar 
masses at the GUT scale are set to zero, mo = 0, and neglecting all the other positive 
but small contributions (proportional to af) we can estimate the minimum mass of 
Mr as 

"R 



Ml ^ -MiiMs) 



al{M, 



X/ 



«i(M3 



(4.6) 



From the pole gluino mass, Mg, we extract the DR mass pO 

M3(M3 



M~g 



1- 



15«3(M3 



A-n 



(4.7) 



For Mg = 200 GeV, we get 1/3(71/3) = 174 GeV" and 



Mr ^ 
Or 



150 GeV. 



(4. 



For larger gluino masses this becomes larger. If we allow for a positive contribution 
from ml then M^ becomes correspondingly larger. From Figs. [l|, |^ we see that this 
value is incompatible with a light sbottom of 0{< lOGeV) in the small tan/5 ~ 10 
region where the above solution is valid. In fact, universality of the squark and 
slepton masses is incompatible with all the values of tan/5. This is because ttiq = 
implies a small stop mass (even smaller than the sbottom one) excluded by the 
current experimental data. 

Let us now consider the case of a light gluino here taken to be 3 GeV. In order 
to obtain chargino and neutralino masses compatible with the experimental data 
we keep the common electroweak gaugino mass M2 = Mi = M1/2 at the GUT 
scale greater than 120 GeV. Then for Mg = 3 GeV and as{mf,) = 0.22 we get 
M^lM^ ~ rrif,) = 2.2 GeV. Evolving this up to the Z-scale using the relation 

as^mb) as^Mz) ' 
we obtain Ms{Mz) = 1.2 GeV which in turn gives from (f4.4|) 

M-,^ ~ mo. (4.10) 

That is compatible with the Mj^ mass of our Fig. |^ for positive rriQ at the GUT 
scale but again is not compatible with the experimental bound on the top squark 
mass (> 120 GeV). Thus we conclude here that a light sbottom mass of order 
~ 0(10 GeV) is incompatible within the MSSM under the assumption of universality 
of scalar masses as well as universality of the electroweak gaugino masses at the GUT 
scale. The gaugino mass universality is not essential. 

4.2 Non-Universal Scalar Masses at Mx 

Analytical solutions of the renormalization group equation ([4.1|) in the case of non- 
universal boundary conditions have been obtained in [^ , under the assumption of a 
small bottom Yukawa coupling and thus small tan f3 values. Even in this approxima- 
tion the results are quite complicated. The term Tr(ym^) of ( [4.3| ) is now non-zero 
at the GUT scale and below. One must thus consider the running of the other soft 
masses as well. The coupled system of differential equations is difficult to solve. 

Since we are interested in solutions of the RGE's even in the large tan /3-regime 
we solve them numerically. Instead of solving the RGE's assuming a specific pattern 
for the soft breaking masses at the GUT scale, we use our results from Figures lb and 
3a for the DR right handed soft bottom mass at the Z-scale and run this up to the 
GUT scale together with all the other masses and couplings. We use two-loop RGE's 
for all the couplings and masses and full treatment of threshold effects [^]. All the 



other parameters have been taken to satisfy the current experimental constraints. 

The dominant effect on the running of Mi is the gluino mass. As we run the 
RGE's up in scale, the gluino mass drives M^ to negative values. The scale where 



M^ becomes negative depends on tan/3 and on the initial M^ . Note that the 
positive bottom Yukawa couphng contribution compensates the negative ones from 
the gluino in ( ^l]) "'. 

M^ < imphes a charge and colour breaking minimum of the scalar potential 
(CCB). Note that the CCB we obtain is above a given scale, i.e. at the weak scale 
or below charge and colour symmetry are restored. Such a CCB at high scales (as 
opposed to present scales [^) is not observationally excluded to our knowledge [pO| 



however it would substantially alter the conventional cosmology at high tempera- 
tures. If colour and charge are broken close to the electroweak scale we would expect 
there to be experimentally observable effects. It is beyond the scope of this letter 
to investigate this in detail. In the following we shall consider charge and colour 
breaking to be excluded for scales Q < I TeV. 

In order to avoid CCB we must go to large initial values of M^ . As we saw 
in our approximate analytical solution, the only possibility then of getting a light 
bottom squark (pole) mass is to go to large values of tan jS. This is confirmed by our 
full numerical calculations. In the following Table we summarize the lower bound 
on tan/5 for a given gluino mass for which there is no CCB at all scales below 
Q < ITeV. In fact there is also an upper bound on tan/5. For large values of tan/? 
the light sbottom squared pole mass becomes negative. In other words the gluino 
contribution to the physical sbottom mass is such that there is always a tachyon. 
The upper bound on tan (3 is also summarized in the Table below. 



M, 


200 GeV 


250 GeV 


300 GeV 


350 GeV 


tan/? > 


19 


23 


28 


33 


tau/S < 


53 


48 


40 


34 



(4.11) 



So for Mg = 200 GeV we only get a light bottom squark without CCB for tan/? > 19 
and we avoid a tachyonic sbottom for tan/5 < 53. For Mg > 350 GeV there is no 
physical light sbottom mass in the MSSM. One could rather expect this result : A 
light sbottom of order of 10 GeV or lighter is potentially unstable under the radiative 
corrections of particles with mass of a few hundred GeV. 

Next, we consider the case where we require that there is no CCB for all values 
Q < Mx, i.e. up to the GUT scale. The corresponding excluded parameter range 
is shown as the large shaded region in Figs, lb and 3a for a heavy and light gluino, 
respectively. For a heavy gluino, we see that in order to retain the conventional 
cosmology we must require tan/5 > 32 and M^ > 90 GeV. For a heavier gluino, 
Mg > 300 GeV^, the entire tan/5 < 50 plane is excluded. For a light gluino (Fig 3a) 
this excluded region is substantially reduced. Also M? (Q) turns negative only very 
close to the GUT scale at about Q ~ 5 x 10^^ GeV thus possibly avoiding most 
cosmological constraints. 



Non-universality affects also the results through the term Tr (Ym) in (4.1) 



In summary, in the case of a heavy gluino in order to have a hght bottom squark 
we typically obtain CCB. The heavier the gluino the lower the scale at which CCB is 
obtained. CCB below 1 TeV is only avoided for large values of tan/?, as summarized 
in ( [4.11| ). In order to avoid CCB altogether we must go to very high values of tan/3 
and Ml as shown in Fig. lb. These constraints are largely avoided for a light gluino, 
as seen in Fig. 3a. 

So far we have not considered radiative electroweak symmetry breaking (RESB) 
[0. This is potentially a very strict requirement. As we saw, for a heavy gluino we 
required large tan/3 in order to obtain a light bottom squark. This leads to a large 
bottom quark Yukawa coupling such that possibly both M^^ and M^ are negative. 
This is inconsistent with electroweak symmetry breaking. A systematic check of this 
constraint is beyond the scope of this letter. However, we would expect it to possibly 
be important. 

5. Light Gluino 

As we have discussed, the introduction of a light gluino naturally allows for a light 
M^ even for small values of tan/?, while avoiding any CCBs. However, it appears 
that such a light gluino is experimentally excluded. In order to discuss this we 
distinguish between a decaying and a non-decaying light gluino. The latter could for 
example be the LSP. We first discuss the case of a decaying light gluino. 

Until fairly recently, there was a window in the search for a decaying light gluino 
after combining several sets of experimental data [^, Q . However, this window has 



now been closed [Q by new data from KTeV |]55[ and from LEP [^ ^|. We do not 
consider it any further. 

Next we consider a stable light gluino. (It would be stable if it were the LSP.) It 
would have been produced in the early universe and would have a non- vanishing relic 
density today. These relic gluinos could bind with nuclei (possibly after forming a 
bound state, such as R^ = gg) leading to anomalously heavy nuclei. The number of 
such nuclei depends on the relic density (which is a function of the self-annihilation 
cross section) and also on the binding potential with nuclei (which depends on the 
scattering cross section with nuclei). The relic density was first considered in [p^j. 



The resulting number of anomalously heavy nuclei present today were shown to be 
excluded by existing searches in [^]. More recently this problem has been revisited 
with more detailed work on the binding potential with nuclei ||3^ , however with the 
same conclusion, excluding a stable gluino [^ . 



In ^U the self-annihilation cross section of the gluinos was reinvestigated. The 
authors concluded that unknown non-perturbative effects could possibly lead to a 
larger cross section and thus a significantly smaller relic density. This could possibly 
avoid the bounds from anomalous heavy nuclei searches. The authors take this as a 
motivation to re-examine bounds from colliders. Using existing analyses from LEP 
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(OPAL 1^) and the Tevatron (CDF |^) they exclude the full range of gluino masses 
homSGeV-lSOGeV.^ 



6. Conclusions 

We have investigated the case of whether a light bottom squark C(< 10) GeV can 
be accommodated in the unconstrained MSSM. In our analysis we have included all 
the relevant one-loop corrections to the physical sbottom pole mass. For tan/3 ~ 15 
these corrections are large and need to be included when determining the physical 
parameters. The main effect is from the heavy gluino mass, but also the trilinear 
coupling Ab leads to significant effects. In this precise framework, we were able to 
extract the running parameters and evolve them up to higher scales with the full 
two-loop RGEs including all threshold effects. In detail we find: 

• If we assume universal scalar masses at the GUT scale (minimal supergravity 
scenario) we find a light sbottom is inconsistent with the experimental bounds 
on the other supersymmetric scalars. Thus in this scenario a light sbottom is 
excluded. 

• A light sbottom can be embedded in the MSSM with non-universal scalar 
boundary conditions at the GUT scale only for specific conditions. For a heavy 
gluino {Mg > 180 GeV) it requires large values of of tan/3 > 30, in order to 
avoid CCB at scales below Q < 1 TeV. This lower bound on tan f3 grows with 
the gluino mass, ( [4.11J ). For each gluino mass there is also an upper bound on 



tan/3 beyond which the light sbottom mass becomes tachyonic, ( 4.11 ). Above 
Mg > 350 GeV a light sbottom is completely excluded. The gluino mass is thus 
restricted to the range 180 GeV < Mg < 350 GeV. 

If we require the absence of CCBs up to the GUT scale the allowed values 
of Mg and tan/3 are significantly more restricted as summarized in Fig. lb. 
Gluino masses above 300 GeV are already excluded. 

A light sbottom could be embedded naturally in the MSSM with a light gluino 
~ 3GeV in a less fine tuned way avoiding also CCB constraints for almost all 
the tan /3 values. However, a light gluino seems to be experimentally unlikely. 



^There is a possible window between 25 GeV — 35 GeV for an "unlikely" set of parameters |4l| , 
which is not of direct interest to our problem. The phenomenological consequences of this window 
have been further explored in |E3. 
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We conclude that a light sbottom hypothesis is not completely excluded in the 
MSSM but it is disfavoured. 



Note Added: After completing this paper, Piep-ph/0008321| by Plehn and Nierste 



was put on the net. It is complementary to our work focusing on effects in the 
B-meson data for specific bottom decays to sbottoms. 
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